home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / sgetri.z / sgetri
Encoding:
Text File  |  1998-10-30  |  2.2 KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSGGGGEEEETTTTRRRRIIII((((3333FFFF))))                                                          SSSSGGGGEEEETTTTRRRRIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SGETRI - compute the inverse of a matrix using the LU factorization
  10.      computed by SGETRF
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE SGETRI( N, A, LDA, IPIV, WORK, LWORK, INFO )
  14.  
  15.          INTEGER        INFO, LDA, LWORK, N
  16.  
  17.          INTEGER        IPIV( * )
  18.  
  19.          REAL           A( LDA, * ), WORK( LWORK )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      SGETRI computes the inverse of a matrix using the LU factorization
  23.      computed by SGETRF.
  24.  
  25.      This method inverts U and then computes inv(A) by solving the system
  26.      inv(A)*L = inv(U) for inv(A).
  27.  
  28.  
  29. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  30.      N       (input) INTEGER
  31.              The order of the matrix A.  N >= 0.
  32.  
  33.      A       (input/output) REAL array, dimension (LDA,N)
  34.              On entry, the factors L and U from the factorization A = P*L*U as
  35.              computed by SGETRF.  On exit, if INFO = 0, the inverse of the
  36.              original matrix A.
  37.  
  38.      LDA     (input) INTEGER
  39.              The leading dimension of the array A.  LDA >= max(1,N).
  40.  
  41.      IPIV    (input) INTEGER array, dimension (N)
  42.              The pivot indices from SGETRF; for 1<=i<=N, row i of the matrix
  43.              was interchanged with row IPIV(i).
  44.  
  45.      WORK    (workspace/output) REAL array, dimension (LWORK)
  46.              On exit, if INFO=0, then WORK(1) returns the optimal LWORK.
  47.  
  48.      LWORK   (input) INTEGER
  49.              The dimension of the array WORK.  LWORK >= max(1,N).  For optimal
  50.              performance LWORK >= N*NB, where NB is the optimal blocksize
  51.              returned by ILAENV.
  52.  
  53.      INFO    (output) INTEGER
  54.              = 0:  successful exit
  55.              < 0:  if INFO = -i, the i-th argument had an illegal value
  56.              > 0:  if INFO = i, U(i,i) is exactly zero; the matrix is singular
  57.              and its inverse could not be computed.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.